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A long-range force acting only between nonbaryonic particles would be associated with a large 
violation of the weak equivalence principle. We explore cosmological consequences of this idea, 
which we label ReBEL (daRk Breaking Equivalence principLe). A high resolution hydrodynamical 
simulation of the distributions of baryons and dark matter confirms our previous findings that 
a ReBEL force of comparable strength to gravity on comoving scales of about lh~^Mpc causes 
voids between the concentrations of large galaxies to be more nearly empty, suppresses accretion 
of intergalactic matter onto galaxies at low redshift, and produces an early generation of dense 
dark matter halos. A preliminary analysis indicates the ReBEL scenario is consistent with the 
one-dimensional power spectrum of the Lyman- Alpha forest and the three-dimensional galaxy auto- 
correlation function. Segregation of baryons and DM in galaxies and systems of galaxies is a strong 
prediction of ReBEL. ReBEL naturally correlates the baryon mass fraction in groups and clusters 
of galaxies with the system mass, in agreement with recent measurements. 

PACS numbers: 98.80.-k, ll.25.-w, 95.35. -|-d, 98.65.Dx 



I. INTRODUCTION 

The ACDM cosmology (the relativistic hot Friedmann- 
Lemaitre model with a cosmological constant and Cold 
Dark Matter, DM) is very successful at matching the 
large-scale distributions of matter and radiation, includ- 
ing the angular Power Spectrum (PS) and polarization 
of the cosmic microwave background and the galaxy cor- 
relation functions and PS, if one allows for mild biasing 
between the mass and galaxy distributions. But there 
are problems with structure formation, and there is a 
possible remedy, ReBEL, a long-range force of attraction 
operating only on the DM. This has been studied in nu- 
merical simulations of structure formation in QtiM ■ Here 
we present new results on structure formation on cosmo- 
logical length scales from joint modeling of baryons and 
DM using the Smoothed-Particle Hydrodynamics (SPH) 
method and the DM TreePM code Gadget2 5], modi- 
fied to take account of the added force on the DM. Our 
two simulations compare evolution of structure with and 
without ReBEL. We also use the Halo Occupation Dis- 
tribution (HOD) framework to compare model predic- 
tions to statistical measures of large-scale structure de- 
rived from recent galaxy surveys. 

We consider ReBEL in the form of an attractive force 
between DM particles alone of the form 
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as in P, 13, Here ar is the physical separation vec- 
tor, where a{t) is the cosmological expansion parameter. 
The screening length, here set to = lh~-'^Mpc (where 
the Hubble parameter is Hq = 100ft- km s~^ Mpc~^), is 
constant in comoving coordinates. The Newtonian grav- 
itational constant is G and the DM particle mass is m. 

Frieman and Gradwohl pointed out the constant /3 that 
measures the strength of the ReBEL force is bounded 
above by its tendency to separate stars from DM 
contrary to the observation that some dwarf spheroidal 
companions of the Milky Way have retained their DM 
halos. Kesden & Kamionkowsky [21 made the excellent 
point that a moderately strong ReBEL, with 13 ^ 0.1, 
separates stars from DM in a satellite galaxy in one-sided 
stream. This is quite contrary to the observed classical 
double stream of stars in the Sagittarius galaxy @ . We 
consider here the stronger ReBEL force (3 = 1. This 
passes the Frieman- Gradwohl test in the cases consid- 
ered in Q . It also passes the Kesden-Kamionkowsky test 
because ReBEL would pull the Sagittarius stars away 
from the DM halo before close passages by the Milky 
Way tidally disrupt the DM-free but still gravitationally 
bound star cluster in the observed two streams. The pro- 
cess is illustrated in Figures 3 and 11 in Keselman et. al. 
Q. Bean et al. :9] show that /3 ~ 1 also so far passes the 
constraints from the cosmological tests. 

Section |TT] reviews issues that motivate analysis of ob- 
servational consequences of ReBEL, with a guide to the 
issues studied in the present paper. The numerical mod- 
elling is described in § IIIII Results are presented in § IIVI 
beginning with a general description of the main results. 
The more elaborate analyses are divided into several sub- 
sections describing the different algorithms used for find- 
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ing halos and voids and the method used to calculate the 
Lya PS. A concludmg discussion is presented in §|Vl 



II. ISSUES AND MOTIVATIONS 

We review here our selection of apparent challenges to 
structure formation that motivate consideration of the 
addition of RcBEL to ACDM. For another overview of 
the situation see [loj . 

The ACDM cosmology tends to over-populate halos 
of galaxies with satellites and voids with dwarf galax- 
ies. The former [ll| - [T5| may be resolved by the effect of 
plasma pressure or winds that leave low mass DM ha- 
los with too few baryons to be observable. ReBEL of- 
fers another stripping effect, the differential acceleration 
of baryons and DM [4], and ReBEL may also promote 
more coniplete merging of small halos [11, H| . The latter 
issue pi3l - [l8j is best expressed as a comparison of counts 
of galaxies actually observed in denser regions to what 
is seen in voids [l^ [13]. In particular, the nearby Lo- 
cal Void contains just two known galaxies [2ll - l23l | , while 
scaling from counts of galaxies observed elsewhere would 
predict the presence of about 15 to 30 galaxies with 
— 18 ^ ^ — 10 in the Local Void. The possible role 
of a ReBEL force in more completely emptying voids is 
explored in p5] and in § IIV CI 

While ACDM correctly predicts the general density 
profile of clusters of galaxies it tends to underestimate 
the concentration parameter, resulting in profiles that 
are too shallow with respect to observations based on 
strong and weak lensing [3, [23 - [26| and X-Ray emission 
[23 ■ We show in S IIV B 31 that ReBEL increases the con- 
centration by earlier assembly of clusters, which is in the 
direction indicated by the observations. 

The ACDM cosmology underpredicts the abundance of 
disk-dominated galaxies [H, [2^ . This is because merg- 
ing accretion of extragalactic debris by DM halos of mass 
10^^ — 10^'^ Mq is predicted to continue to modest red- 
shifts, causing thin disks that formed earlier to thicken 
into bulges by gravitational heating [30l - l33| . The pre- 
dicted merging also is difficult to reconcile with the in- 
sensitivity to environment of the correlations among lu- 
minosity, radius, velocity dispersion and color in early- 
type and late- type galaxies [34l-l37| . We show in § IIV B 21 
that ReBEL promotes significantly earlier assembly of 
galaxies, which is in the direction wanted to relieve these 
problems. 

The earlier structure formation under ReBEL would 
promote earlier reionization. Polarization measurements 
of the cosmic microwave background radiation indicate 
reionization of hydrogen by redshift 11, implying early 
structure formation [38, 39]. Galaxy/star formation in 
the ACDM with as ~ 0.9 could be early enough to ac- 
count for such early reionization. However, extreme ef- 
ficiency of UV photon production by early structures is 
needed for erg ~ 0.8, the value favored by WMAP data 
[ioj . ReBEL may relieve this condition. 



Earlier merging would lead to stronger clustering of 
more massive galaxies now and the clustering of the in- 
tergalactic medium at redshift z ^ 3 observed as the 
Lyman- Alpha forest. The galaxy two-point correlation 
function and PS found in [^ in and here in IIV A 31 seem 
to be consistent with what is observed within the usual 
uncertainties of biasing. Dark matter halos are more 
massive in ReBEL (§ IIVB 1[) . and the difference from 
the standard model is more pronounced at redshift z = 1 
than the present. Remaining to be explored is the ef- 
fect on the galaxy luminosity function. We conclude in 
§ IIV A 21 that the model parameters for the intergalac- 
tic medium at x ~ 3 can be adjusted within the avail- 
able constraints to fit the one-dimensional Lya PS about 
equally well in the standard model and in ReBEL with 

Groups of galaxies in which the virial temperature is 
large enough that intragroup plasma would be expected 
to be detectable in X-ray emission are observed to con- 
tain a significantly smaller baryon mass fractions than 
the cosmic mean (e.g. [4T| - |45j ). This may be a result 
of the tendency of ReBEL to segregate DM from bary- 
onic matter, as shown in [46l and further analyzed in 

§IIvHl 

The more rapid structure formation that seems to be 
indicated by the phenomenology we have reviewed could 
follow from at least three adjustments of ACDM: a de- 
parture from near scale-invariant adiabatic Gaussian ini- 
tial conditions, a modified gravity with a departure of 
the inverse square law that preserves the weak equiv- 
alence principle [1^, or, in general relativity theory, a 
long-range force that operates only on the DM parti- 
cles. We consider this last option. Superstrin^s offer 
scenarios for such a force in the dark sector [49l - l51| . The 
force would cause the accelerations of cosmic objects to 
depend on their relative DM content, in what may be 
termed a breaking of the weak equivalence principle, or 
ReBEL (daRk Breaking Equivalence principLe). Bovy 
and Farrar 52] point out that if the DM interacted with 
visible matter strongly enough to allow laboratory DM 
detection then ReBEL could induced a fifth force in the 
visible sector that is much stronger than allowed by the 
Eotvos laboratory constraint. It should be noted, how- 
ever, that a laboratory DM detection need not vitiate 
ReBEL, for there may be several types of DM in addi- 
tion to massive neutrinos, including one with a significant 
interaction with the visible sector but not ReBEL , and 
another that is coupled to ReBEL but not the visible 
sector [iiil. 



III. NUMERICAL MODELLING 

A. The simulations 

We have modified the Gadget2 code to include, in 
addition to gravity, the force in Eq. [T] operating be- 
tween the DM particles only. The simulations include 
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baryons, (3=1 baryons, (3=0 DM, P=l DM, P=0 




X [h 'Mpc] X [h 'Mpc] X [h 'Mpc] X [h 'Mpc] 



FIG. 1: Particle distributions in a slice 2h ^Mpc thick. Only a random subsample, containing 10% of the particles is shown. 



gas (baryonic) particles subject to pressure and gravita- 
tional forces but not the ReBEL force. The hydrody- 
namic equations are solved using the particle-based SPH 
method in Gadget2. The modification is such that the 
error in the calculation of the total force is the same as in 
the standard, public Gadget2. Appendix El describes the 
details of the modification. The modified code is tested 
as described in appendix IB] 

We obtain for comparison two simulations, with and 
without ReBEL. The cosmological background is the 
same in the two simulations, with density parameters 
^DM = 0.23, riBAR = 0.046, riA = 0.72, with Hubble pa- 
rameter h — 0.7, all consistent with [55;]. Encouraged by 
our previous results T, "3] we adopt the screening length 
rs = 1 h~^Mpc and RcBEL force parameter /? = 1 in Eq. 

m 

The simulation box is GOh^^Mpc on the side. It con- 
tains 128'^ DM particles and an equal number of baryonic 
particles. The DM particles have mass 6.57 x 10^ h^^M©, 
the baryon particles 1.31 x 10^h~^MQ. The bayonic 
equation of state is that of a perfect gas with 7 — 5/3, 
neglecting energy sinks and sources. The specific entropy 
is conserved except at shocks. The initial entropy corre- 



sponds to a temperature of 5 K at the starting redshift 
of the simulation, z = 1000. Shock waves generated in 
collapsing structures heats the gas to the much higher 
temperatures noted in the next section. 



Our two simulations start from identical initial condi- 
tions at z = 1000, when the ReBEL force is assumed to 
start operating [511 ] . The initial conditions a re g enerated 
from the linear PS of the standard ACDM [5g] cosmol- 
ogy for our choice of parameters. We use the Zel'dovich 
displacements of points and their peculiar velocities in 
an initially uniform body-centered cubic lattice. The ini- 
tial density fiuctuations are normalized so the linearly 
extrapolated (without ReBEL ) rms value of the density 
fiuctuations at z = in spheres of radius 8h~^Mpc is 
erg = 0.8. 



To identify DM halos in the simulation we apply the 
"Friends-of- Friends" (EOF) algorithm [13]. The inter- 
particle linking length is b — 0.2 in units of the mean 
particle separation. 
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B. Galaxies and Halos in the simulations 

To compare results from the simulations to the ob- 
served large-scale distribution of galaxies we need a recipe 
for identifying mock galaxies. The most thorough way is 
to trace halos in the simulation forward in time and as- 
sign them galaxies according to semi-analytic galaxy for- 
mation recipes (e.g. [SS]), but this requires simulations 
with larger resolution than presented here. We use in- 
stead the s imp ler Halo Occupation Distribution (HOD) 
framework [59| , in which a DM halo with mass Afj^ above 
a threshold M^j^ contains a central galaxy, and, above a 
larger threshold Mi, satellite galaxies. 

We adopt HOD parameters such that galaxies in the 
simulations represent the observed galaxies more lumi- 
nous than Li,/2.b, corresponding to absolute magni- 
tude Mr = —19.5, with the observed number density 
O.Olh^Mpc"^ adopted from table 2 in [gO]- In the con- 
ventional case, /? = 0, the HOD parameters also produce 
a satisfactory fit to the observed number density of voids 
above a mass threshold of 10"'^° h~-'^M0 and the distri- 
butions of galaxies within them [l^]. The /3 = 1 case 
requires a modification of the HOD prescription to fit 
the voids, as we will describe. 

In the P — Q case the number of satellite galaxies as- 
signed to a halo of mass is drawn from a Poison 
distribution with mean 




The three parameters A/^.^^, Afi, and a are fixed to 
match the galaxy number density and correlation func- 
tion. Each satellite is placed on a randomly chosen halo 
DM particle. 

To match both the the properties of voids and the dis- 
tribution of bright galaxies in the /3 = 1 case we must 
modify the HOD prescription. We adopt a broken power- 
law model for the number of satellites in a halo, 

f fxr-V if M <M , 

i (mt) + (tT^J otherwise, 

which is continuous at the mass scale M^. Also, we need 
a more extended distribution of satellite galaxies than 
results from the jS — Q prescription. We achieve this by 
increasing the satellite distances from the halo center by 
a constant multiplicative factor, C, relative to the DM 
distribution in the host halo. This produces a satellite 
distribution similar to that of the /3 = case jFIl, !62| . 

IV. RESULTS 

Fig. [T] shows a representative distribution of the gas 
and DM particles in a slice 2h~^Mpc thick through the 
simulations at three redshifts, as indicated in the figure. 
Comparison between the DM distributions (two columns 



to the right) with ReBEL (/3 = 1) and the standard 
model {P = 0) clearly shows the ReBEL enhancement 
of small-scale clustering of the DM: voids are emptier 
and the presence of individual chimps more pronounced. 
There also is a stronger redshift evolution in the DM 
particle distribution for P — 1. (Here and throughout 
we use coordinates comoving with the general expan- 
sion.) In the /3 = simulation the baryonic component 
closely follows the DM on scales larger than the Jeans 
length, as expected. Also as expected is the behavior 
for /3 = 1, where the baryons are affected by the ReBEL 
force only indirectly, through the gravitational attraction 
of the more strongly clustered DM. The differences be- 
tween the baryon distributions in the f3 = 1 and /3 = 
simulations in the two columns to the left, though im- 
portant, are indeed less pronounced than in the DM dis- 
tributions. 
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FIG. 2: Scatter plot of temperature versus density of particles 
in the 13 — (left panels) and /3 = 1 (right) simulations, at 
redshifts z — (top panels) and z — 3 (bottom) 

Fig. [2] shows a scatter plot of temperatures and over- 
densites of 10^ SPH gas particles in the two simulations 
at two redshifts. The top-right corner of this T — p plane 
is more densely populated for /? = 1, and the difference 
is more pronounced at z = 3: in ReBEL more baryons 
are in dense hot regions. Apart from that the diagrams 
for ACDM and ReBEL are similar, consistent with the 
rough similarity of the baryon distributions in Fig. [T] 



A. Correlation Functions and Power Spectra 

1. The distributions of baryons and dark matter 

The clustering properties of the particle distributions 
measured by their PS are shown in in Fig. [3] at several 
redshifts. Thin lines correspond to the standard scenario, 
/? = 0, while the thick lines show the effect of ReBEL with 
13 = 1. As expected, the PS in the two simulations are 
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FIG. 3: PS of the baryon and DM distributions for /3 = and 
/? = 1 at the redshifts indicated in the figure. 



almost the same on scales large compared to the screen- 
ing length, = lh~^Mpc, while ReBEL considerably 
increases the clustering of the DM on smaller scales. In 
the standard scenario the baryon PS (thin dashed) closely 
follows that of the DM (thin solid) at all redshifts, while 
in the RcBEL scenario the baryon PS (thick dashed) falls 
below the DM curve (thick solid). The baryons in the 
ReBEL simulation are significantly more clustered than 
in the standard model at higher redshifts, but at z = 
the baryons in the standard and ReBEL scenarios have 
relaxed to quite similar PS. This evolution is illustrated 
in another way by the autocorrelation functions in Fig. 
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Galaxy redshift surveys provide redshift distances that 



loair/Mpc/h 1 

FIG. 4: Autocorrelation functions for the two components of 
the two simulations at the current epoch and redshift z = \. 



differ from real distances by the radial velocity term, 
breaking the isotropy of the position correlation function, 
thus offering a probe of the peculiar velocity field and 
through that the cosmological model. Fig. [S] shows the 
redshift space correlation function f (^p, tt) as a function 
of the redshift space separation perpendicular to the 
line of sight and tt parallel to the line of sight [63^, '63| , but 
here computed for the DM in the standard and ReBEL 
cases. Random motions on small scales cause the spike 
at Tp « 0. It is more pronounced for /? = 1, reflect- 
ing the stronger small-scale clustering. On scales larger 
than r^= \ h^^Mpc, the contours have a similar pattern. 
We cannot explore the anisotropy of the redshift space 
correlation functions of halos and mock HOD galaxies 
because the anisotropy is sensitive to peculiar motions 
and our simulations are of insufficient dynamical range 
to model these motions properly. 



2. The Lya jorest 

The Lya forest in QSO spectra is caused by resonant 
scattering of QSO light by diffuse atomic hydrogen along 
the line of sight. This is an important tracer of the dis- 
tribution of intergalactic baryons at high redshift f65|. 
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FIG. 5: Contours of the DM redshift space correlation func- 
tion ^"(rpjTr) as function of perpendicular (rp) and parallel 
(tt) separations. The contours are spaced by intervals of 0.2 
in log(^''), starting with log(^^) = —1.9 at the outermost con- 
tour. 



To estimate Pgas as a function of position in the sim- 
ulations we interpolate the baryon particle positions by 
a cell-in-cloud (CIC) procedure to derive the gas den- 
sity field on a 256'^ uniform cubic grid in the simulation 
box. This grid has a Nyquist frequency close to the up- 
per limit beyond which the flux PS is contaminated by 
metal lines. The density field is then deconvolved by the 
CIC kernel sine (kxL/Ng) sine (kyL/Ng) sine (kzL/Ng) 
where L is the box length, and Ng is the number of cells 
on the side. The density field along 256^ sight lines is 
then smoothed with a ID Gaussian window of width aj. 
This smoothing is temperature-dependent and is intro- 
duced in order to mimic temperature broadening [6^ . It 
will be treated here as a free parameter. The resolution 
of the simulations is insufficient to model the PS properly 
at wavenumbers close to the upper limit in EqlHl so we 
degrade the observations to conform with the resolution 
of our simulations. We convolve the observed PS with a 
Gaussian window [69], 



Ppik) =^ Ppik) e^p{k^/kl 



(7) 



The standard ACDM scenario matches well the observed 
one-dimensional (ID) PS derived from the forest [gl]. We 
present here a preliminary assessment of whether or not 
the ReBEL scenario can also match the ID PS. This is an 
interesting test because it probes the baryon distribution 
at epochs where the ReBEL scenario significantly differs 
from the standard model (Fig. . 

The normalized transmitted fiux (hereafter fiux) in the 
Lya absorption spectrum is 



F — exp(— t) 



(4) 



The optical depth r depends on the mass density pgas 
in baryons, plasma plus neutral, where a is a parameter 
in the range 1.5 to 2, depending on the photoionization 
history, and the proportionality factor A, which depends 
on the intensity of the ionizing background [13] , is given 
by M 



Aiz) 



0.61 



(5) 
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0.02 
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1.5 X 104 K 



Here J7bar is the baryon density parameter and Fphot is 
the photoionization rate per hydrogen atom. 

We consider the Lya ID PS of F for wavenumbers in 
the range 



0.3h Mpc" 



fc < 4h Mpc^ 



(6) 



The lower bound is set by fiuctuations in the continuum 
69] and the upper bound by contamination by metal lines 
20, El- We compare the simulations to the observed Lya 

ID PS at redshift z « 3. 



where ka- = VS In2/FWHM. We use a FWHM equiv- 
alent to the mean separation between gas particles. In 
addition, to avoid dealing with continuum fitting effects, 
we truncate the transmitted flux 66i at F = 0.8. 




log[k/Mpc 7h] 



FIG. 6: The ID Lya flux power spectra computed from the 
/? = and P — 1 simulations are shown as the dashed and 
solid lines respectively. The observed PS is shown as circles 
overlaid with crosses representing the error bars and bin sizes. 
The vertical dotted lines correspond to the bounds in Eq. [B] 

The dimensionless measure k Ppik) of the observed ID 
PS of the Lya transmitted flux compiled in [t^ is plot- 
ted as the open circles in Fig. |6l The curves show the 
ID PS computed from the two simulations. The param- 
eters in these simulated PS are, for /3 = 0, a = 1.99, 
A = 0.48, and aj = 0.32 fi-^Mpc, and, for /3 1, 
a = 1.66, A = 0.71, and cr/ = 0.35h~iMpc. For mean 
ICM temperature T — 1.5 x 10'' K these parameters give 



photoionization rates Fphot = 6.9 x 10 s for 
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and 4.6x 10 ^^s ^ for /3 = 1. The values of cr/ correspond 
to 18 h Mpc~^ for the broadening parameter used 

in [66|, shghtly larger than the maximum value adopted 
in that work. 

Fig. |6] shows that both ACDM and ReBEL can agree 
with the observations within the limits of the present 
analysis. We caution, however, that this is a prelimi- 
nary test because our simulations neither include a direct 
treatment of photoheating and photoionization nor they 
are of sufficient resolution to resolve the Jeans mass. 



3. Clustering of galaxies 

Galaxies form in DM halos where the gas density can 
reach high enough values to cool, contract, and form 
stars. Therefore, in contrast to the Lya forest, which 
traces the baryons, the distribution of galaxies is dictated 
by the DM density field. 



4.5 




-1.5 -1 -0.5 0.5 
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FIG. 8: Galaxy PS from HOD galaxies in the simulations and 
in the SDSS galaxy survey. 
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FIG. 7: Correlation functions for the positions of the HOD- 
placed galaxies in the P — and /? = 1 simulations are shown 
as the thin dashed and thin solid lines respectively. The ob- 
served correlation function of galaxies is derived from the 2dF 
and SDSS surveys. 

In the simulations we populate halos with "galaxies" 
according to the HOD framework described in § IIIIBI 
The HOD parameters for the (3 = simulation are 
a = 1.13, M„,„ = 3 X 10", and M, = 8 x lO^^h-iMQ. 
They are similar to the ones derived in (60j . The parame- 
ters for the /3 = 1 simulation are 1.3 x 10^^ h-^Mg, 
a = 1.13 (as for /3 = 0), 7 = 1.8, M^.„ = 2.2 x lO^^, 
Afj = 6.3 X 10^2 jj-i^^^ ^jjg p =Tsimulation the 
multiplicative factor, C.^ by which satellite distances from 
the halo center are increased relative to the DM distri- 
bution in the host halo is C = 4. 

The solid and dashed curves in Fig. [7] show the corre- 
lation functions of these galaxies in the two simulations. 
The thick solid line is the fit ^ = (r/5.1 h-iMpc)-^-^ to 
the observed correlation function extracted from Fig. 11 



in [59j . In the range of separations shown here the corre- 
lation functions in the two simulations, with the adopted 
HOD parameters, agree reasonably well with each other 
with what is observed. The alternative view from the 
galaxy PS is shown in Fig. HI The observed galaxy PS is 
taken from the SDSS galaxy redshift survey (table 2 of 

An important test not examined here is the shape 
of the correlation function of the simulated galaxies on 
smaller scales, where ReBEL may be expected to have a 
larger effect. That will require simulations with better 
resolution. 

In the range of galaxy separations we can explore the 
fit of the position correlation function of the ReBEL sim- 
ulation HOD galaxies to the observed galaxy correlation 
requires HOD parameter C = 4. That is, ReBEL requires 
that the distribution of satellite galaxies is considerably 
more extended than the DM in the host halo. This is in 
contrast to the standard cosmology in which only a mild 
bias, or none at all (as in this work), is needed to match 
the galaxy correlation function. In high resolution sim- 
ulations of standard ACDM the distribution of the DM 
sub-halo population is more extended than the distribu- 
tion of diffuse DM in the host halo That is, HOD 
satellite galaxies which trace the general DM distribution 
cannot be associated directly with the sub-halo popula- 
tion. We may expect that sub-halos in the ReBEL sce- 
nario are also more extended than the DM distribution 
in the host halo, but that remains to be checked, along 
with the possibility that satellites galaxies could be more 
directly related to sub-halos in the ReBEL scenario. 
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FIG. 9: Mass function of halos in the /3 — and 13 = 1 sim- 
ulations at redshifts z = and z = 1, for DM particles only. 
The vertical dashed line marks a mass of a halo containing 20 
DM particles as identified by the FOF algorithm. 



about two at M ~ lO^^h-^Mo, roughly the DM mass 
of the Milky Way [13] ■ Translating this to the predicted 
effect on the galaxy luminosity function does not seem 
likely to be reliable at the present state of the art. 

The difference between the mass functions is greater 
at redshift z — 1, and the mass function expressed in 
comoving coordinates evolves more slowly at z < 1 in 
the ReBEL scenario. As we discuss next, this agrees 
with the earlier completion of merging in ReBEL. 

li ■ . ■ . . 1 




12 12.5 13 13,5 14 14.5 
logl^^final'^O*"'] 

FIG. 10: Mean fractions of mass assembled by redshifts Z = 
1 and Z = 3 as functions of the present halo mass in the 
simulations with and without ReBEL. 



2. Halo Merging History 



B. Properties of halos 

1. The halo mass function 

To explore the abundance of DM halos as a function of 
their mass we select FOF groups containing more than 
20 DM particles each, giving a minimum halo mass of 
1.3 X lO^^h~^M0. The mass functions expressed as the 
number density of halos per logarithmic mass bin are 
plotted in Fig. [U] The la error bars are estimated by the 
Bootstrap resampling procedure described in l79j . 

The mass functions in the two simulations are nearly 
the same at AI sa 10^^ h^^M©, which is close to the min- 
imum halo mass we can identify. The difference between 
mass functions also is small at M sa lO^''h^^M0, about 
the mass at 2 = at which the virial radius is equal to 
the screening length = 1 h~^Mpc. At the still larger 
masses characteristic of clusters of galaxies the ReBEL 
force is less effective and we may expect the mass func- 
tion to be less sensitive to /3. Our computation box size 
is too small to explore this. 

The mass functions in the ReBEL and standard ACDM 
simulations at the present epoch differ by a factor of 



A present-day halo went through a time of intense 
merging activity when it acquired a significant fraction 
of its present mass. For a measure of when this hap- 
pened wc identify in each DM halo at z = the particles 
belonging to the most massive progenitor (MMP) at red- 
shifts z — 1 and 2 = 3. Fig. [TU] shows the average ratios 
of the mass of the MMP to the present halo mass, M^.^^^j, 
as a function of M,. , . The la error-bars are estimated 
by Bootstrap re-sampling. 

The figure indicates that 2: = 1 a halo with the present 
mass lO^^h^^M0 typical of the Milky Way has assem- 
bled on average 75% of its mass in the ReBEL model 
with l3 — 1, and about half that in the standard model 
with /3 = 0. At2 = 3a halo now typical of the Milky 
Way has assembled 40% of its mass in the ReBEL model, 
while at /? = the MMP typically contains fewer than 
the limiting 20 particles needed for the FOF halo identi- 
fications. Similar results are reported by [8l|. This effect 
is strongest for halos with masses typical of galaxies and 
it nearly disappears at masses typical of rich clusters, 
where it will be recalled that the virial radius exceeds 
the screening length . As discussed in § |TT1 the earlier 
assembly of DM halos in ReBEL with /3 = 1 may play an 
important role in resolving issues of galaxy formation. 
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3. Halo Mass Profiles 

Here we consider the DM mass density as a function 
of radius in halos at z = that satisfy the convergence 
criteria [82j that the two-body relaxation time within the 
studied radius is much longer than the Hubble time, the 
smoothing length of the particles does not allow two- 
body accelerations to be larger than the mean field ac- 
celeration, and there are at least 50 particles within the 
innermost radius studied. That leads us to examine ha- 
los containing each at least 2500 particles identified in the 
FOF algorithm. For each halo we calculate the spheri- 
cally averaged mass profile centered on the most bound 
particle in the halo. Here and in what follows we present 
the mean density p{< r) of mass within radius r, rather 
than the density at radius r, to suppress the noise from 
the limited numbers of particles. 



M=2.7x10^'h"''Mj^ 

op=i ^bi, 
xp=0 ^ 
-NFW fit ^ 


M= 8.5x1 O^V^Mq 


14x10^^h"^Mg 


M=26x10^^h"^M|_, 



-1 -0.5 



-1 -0.5 

log[r«„J 



parameter 

Fit Error ex log 



^ J (logp/PNFw)^ dr 



(9) 



The integration is from rmin to r^^g , p is the measured 
density within r, p^p„ is the fitted NFW density, and 
AR = r^gg — fmin- The second measure is a 'Relaxed' 
parameter. 



Relaxed = log r^gg / r^^ 



(10) 



where r^^ is the distance of the center of mass of the FOF- 
group from its most bound particle. The more spherical 
the halo the larger its 'Relaxed' parameter. 
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FIG. 12: Correlations of halo properties for standard and 
ReBEL model halos. 



FIG. 11: Mass density runs at 2 = in both simulations 
for four randomly chosen halos at the indicated masses. The 
short vertical lines mark the innermost radii that could be 
studied within the convergence requirements. The curves are 
fits to the NFW form. 

Mean density profiles measured for four halos in each 
simulation are shown in Fig. 1111 The inner vertical lines 
attached to each profile indicate the minimum radius, 
rmin, above which the convergence criteria are satisfied. 

The curves in Fig. [11] are integrals of the NFW form 

M 



For each halo in the figure we compute the radius r^^^ 
within which the mean halo density is 200 times the crit- 
ical density, p^ . Then the concentration parameter C^p^ 
is adjusted to fit the mass density run (where p^ is re- 
quired to make the mean density within r^^^ equal to 
200pJ. 

We use two measures of the performance of the NFW 
profile under a ReBEL force. The first is a 'Fit Error' 



Fig. [12] assesses the correlations between pairs of these 
halo structure parameters at z = 0. Circles and crosses 
correspond to /3 = 1 and (3 = 0. The three panels in the 
bottom row are scatter plots of the 'Fit-Error' parame- 
ter versus the halo mass M, the concentration parameter 
Cnfw Eq. m and the 'Relaxed' parameter in Eq. [101 
These panels indicate the NFW profile about equally well 
matches the halo mass profile in both simulations. This 
agrees with the findings of [8I1] , who used a different mea- 
sure than Eq. [S] These scatter plots also make the pos- 
sibly important point that the ReBEL scenario produces 
halos that are more relaxed and more concentrated. 



4- Baryonic mass fraction m halos 

The baryon mass fraction is the ratio of the mass 
density in baryons to the mass density in baryons plus 
dark matter, normalized to unity at the cosmic mean 
ilb/i^m- We consider here the measured [H,!!^ and com- 
puted values of /& in galaxies and systems of galaxies. 

In the largest clusters of galaxies fb is close to unity 
0, HE]; these systems are tightly enough bound that 
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FIG. 13: Baryonic mass fraction versus DM halo mass. Top 
and middle panels, respectively, correspond to masses within 
r^oo and r^gg. Solid curves in the top panel are ilcr limits on 
the mass fraction from recent observations. The bottom panel 
shows the mean gas temperature within r^gg versus the DM 
halo mass within r^gg ; the solid line is the virial temperature 
relation. 



they have resisted large loss of baryons. The heavy blue 
solid curves in the top panel in Fig. [T3] are the ±lcr 
Hmits on observed baryon fractions in lower mass halos. 
This includes the mass in X-ray emitting plasma^and the 
mass in stars, as given in parametric form in 186']. Two 
adjustments are made to these components. First, the 
stellar mass is increased by 11% to take account of the 
contribution from intracluster light [87-89]. Second, ff, is 
reduced by 10%, as a method to overcome the systematic 
errors arising in our simulations, which do not model 



radiative losses and star formation processes, which could 
increase the mass fraction in the simulated haloes by 10% 
[qoI . [qH . Baryonic energy dissipation through radiation 
increases the total baryonic fraction in haloes because gas 
flows deeper into haloes to replace cooled gas. 

The top and middle panels in Fig. [13] show the com- 
puted means of the baryon mass fractions /{, within dis- 
tances r^gg and r^gg from the most bound particle in halos 
at z = in our simulations with /? = (black crosses) and 
/3 = 1 (red circles). The horizontal axis is the DM halo 
mass within r^gg for the top panel and r^gg for the center 
panel. The vertical dotted lines mark the mass in a halo 
containing a minimum mass equal to the sum of 500 DM 
and 500 baryonic particles, the threshold for reliable esti- 
mation of mass fractions in simulations suggested by the 
convergence studies of [Q^, ■ The computed /& in our 
standard ACDM simulation agrees with previous studies 
[q^ . [q^ I . We find that the variation of ft with redshift 
is weak in both the (3 = 1 and (3 = simulations (c.f. 

The bottom panel of the figure shows the mean gas 
temperature within r^gg as a function of the DM mass 
within r^gg . The solid line is computed from the virial 
relation {T ^ GM/r^gg). Halos in the ReBEL (circles) 
and standard (crosses) simulations both closely follow the 
virial relation. 

The top panel of Fig. [Ol illustrates the well-established 
evidence that in galaxies and groups of galaxies roughly 
half the baryons are missing. At halo mass ~ IO^^Mq 
the plasma temperature is '-^ 10^ K, cool enough to es- 
cape detection as an X-ray source. In galaxies at this 
mass and lower the missing baryons could be in the DM 
halo as a galactic corona j97| . At larger mass the missing 
baryons have to have been removed from their DM ha- 
los, perhaps by galactic su perwind s driven by supernovae 
and active galactic nuclei [98l-[l0H. The ReBEL scenario 



adds another possibility, that differential acceleration of 
baryons and DM has separated some baryons from their 
halos. 

The weight of evidence for or against ReBEL from the 
baryon mass fraction /{, in galaxies and systems of galax- 
ies depends on the quite uncertain astrophysics of baryon 
winds. If these processes were effective in suppressing 
then ReBEL could make /h unacceptably small. If winds 
proved to be ineffective it would argue for ReBEL at 
about the parameter values for the DM force we have 
adopted, for the observed and computed baryon frac- 
tions agree in our ReBEL simulation. The ReBEL pre- 
diction for clusters, which have close to the cosmic baryon 
fraction, requires larger simulations. The blue filled cir- 
cle in the middle panel in Fig. [13] shows that the ob- 
served baryon fraction in the Milky Way is well below 
the ACDM case and agrees with ReBEL. This is partic- 
ularly interesting because the relatively low mass of the 
central black hole in our galaxy might argue against re- 
moval of the baryons by a superwind. But the Milky 
Way is at about the mass threshold where the missing 
baryons may be in a corona cool enough to be detected 
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as an X-ray source. 

BAR 



TABLE I: Counts of Voids 
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H = log[l+8] ^l = log[l+8] 



FIG. 14: Probability distributions of the densities in baryons 
and DM smoothed by a top-hat window of radius 1.5 h~^Mpc 



C. Void properties 

Fig. [14] shows the effect of ReBEL on the probability 
distribution functions of the mass densities in baryons 
and DM averaged within randomly placed top-hat win- 
dows of radius 1.5 h~^Mpc. The densities are normalized 
to the cosmic means at the computed redshifts, z ~ 1 
and z = 0, as the density contrast S — p/ p ~ I. ReBEL 
more efficiently evacuates low density regions, as previ- 
ously demonstrated {l|. Here we explore in more detail 
the effects of ReBEL on voids. 

We define a void in the simulation as a spherical region 
of radius R with average density contrast 6v We iden- 
tify voids in the distribution of particles in a simulation 
as follows. We derive a density field by CIC interpola- 
tion of the points to a cubic grid with 128 cells on the 
side. The density field is then smoothed with a Top-Hat 
window with radius large enough that there are no grid 
points with density contrast less than 6v The width of 
the Top-Hat window is then gradually decreased until a 
single grid point with 6 ^ Sv is encountered. This grid 
point is taken to be the center of the largest void and 
its size is the window radius R. The filtering radius is 
further decreased to locate successively smaller voids, ex- 
cluding thos e ove rlapping with the larger, already iden- 
tified, voids (lO^ . 

Table U lists counts of voids in the ACDM and ReBEL 
simulations for two ranges of values of the void radii 
and the threshold density contrast Sv ReBEL consid- 
erably increases the numbers of voids, as also illustrated 
in Fig. [Ml and the effect on the number of larger voids is 
even more pronounced at the lower threshold = —0.9. 
The number of smaller voids is smaller at = —0.8 be- 
cause many of the small voids at = —0.9 become parts 
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FIG. 15: Void mean density within radius r at redshift z = 
for large (bottom) and small voids (top). The void radius is 
R. 



of larger voids at — —0.8. 

We represent the mass distribution within a void as 
the density contrast computed from the mass contained 
within distance r from its center. This measure is less 
noisy than the local mass density profile. Fig. [TS] shows 
means of the void profiles for small (top panel) and large 
voids (bottom). ReBEL reduces the DM density in the 
central regions and produces a steeper density gradient 
at the void edge. Both effects are more pronounced in 
smaller voids, consistent with [l03l |. This could ease the 
tension between the observations that indicate that small 
void profiles have a sharper transition into the outer re- 
gions than what appears in ACDM simulations (l02j . 

The DM particle mass in our simulations is too large to 
allow an analysis of the distribution of low mass halos in 
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voids. We can note, however, that the distinct suppres- 
sion of the mean DM density in voids is in the wanted 
direction to suppress the numbers of void galaxies, while 
the lesser suppression of baryons may be more relevant 
for the density of diffuse plasma in voids. 




loglr/RJ 



FIG. 16: Mean density contrast of HOD galaxies in voids. 

Figs, [in] and [T7] show checks that ReBEL can match 
the void size distribution and the distributions of bright 
galaxies in voids as well as the standard model. The 
galaxies are modeled for absolute magnitude limit —19.5 
in the r band. The HOD parameters were tuned so that 
similar void properties for these galaxies are obtained in 
RcBEL and the standard cosmology. The vertical axis 
in Figs, [in] is the mean HOD galaxy number averaged 
within radius r and scaled to the void radius R. Fig ll7l 
shows the number of voids in our simulations as a func- 
tion of the void radius, using 5y = —0.8. We conclude 
from the consistency of results for /? = and /? = 1 in the 
statistics in Figs. [161 and [TTl that at the level of our simu- 
lations the bright galaxy distribution does not a usefully 
discriminate between ACDM and ReBEL. 



V. DISCUSSION 

We have considered the effect of ReBEL, a long-range 
force of attraction acting on the dark matter alone. The 
force is modeled by Eq. [I] with (3 — 1, meaning the 
strength is the same as gravity, and screening length 

= 1 h~"'^Mpc, meaning it has no influence on the cosmo- 
logical tests but can have interesting effects on galaxies. 

ReBEL with the parameters adopted here satisfies the 
known constraints from properties of satellites of the 
Milky Way, as discussed elsewhere [2i, Q, but this 
certainly requires closer analysis. Within the frame- 
work of the Halo Occupation Distribution (HOD) model 
for galaxies we have shown that ReBEL fits the mea- 
sured galaxy correlation function at separations greater 
than about one megaparsec as well as does the standard 
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FIG. 17: Number of voids smaller than R and larger than 
4h"^Mpc. 

ACDM cosmology. Checking galaxy clustering on smaller 
scales will require larger simulations. 

We have also concluded that ReBEL can account for 
the properties of the Lyman-Alpha forest about as well 
as ACDM, at the level of semi- analytic methods [66lll04l |. 
Again, this requires deeper analysis, including resolution 
of the Jeans mass of the intergalactic medium and an 
explicit treatment of photoionization and photoheating 
by the ionizing radiation background. 

Fig. [131 shows that ReBEL offers a natural explanation 
of the missing baryons in galaxies and groups of galax- 
ies. This may be important because the hypothesis that 
superwinds removed the baryons from isolated galaxies 
and loose groups does not have a secure theoretical or 
empirical basis. 

Fig. [TO] shows that ReBEL substantially advances the 
redshift of assembly of the bulk of the mass of a galaxy. 
This may be important in resolving the puzzle of disk- 
dominated galaxies |l05| . In ACDM inefficient star for- 
mation allows development of dwarf pure disk galaxies 
[l06!|. but that approach would make the disk of the 
Milky Way unacceptably young [107]. ReBEL would help 
by assembling the parts of a galaxy earlier, and assem- 
bling a lower baryon mass fraction. Much better res- 
olution simulations will be needed to explore whether 
ReBEL can promote formation of significant numbers 
of bulge-free spiral galaxies with acceptable galaxy ro- 
tation curves, along with appreciable numbers of large 
red galaxies at high redshifts. 

Earlier structure formation may also be relevant to the 
measurements of polarization of the cosmic microwave 
background radiation that indicate a significant reioniza- 
tion by redshift z ~ 11, implying early structure forma- 
tion. Star formation in the ACDM with cts ~ 0.9 could 
be early enough to account for such early reionization, 
but extreme efficiency of UV photon production by early 
structures is needed for cg k. 0.8, the value favored by 
WMAP data [i^l . ReBEL may relieve this condition. 
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Earlier assembly is accompanied by more massive DM 
halos at the present epoch. The effect is not large, how- 
ever, as illustrated in Fig. [9l and it would seem to be a 
difficult challenge to translate this into a test of ReBEL 
from the galaxy luminosity function. 

Figure [15] shows the depression of the DM mass den- 
sity within voids, and the lesser suppression of the baryon 
mass density. This is in the wanted direction of resolving 
the void over-population problem, but larger simulations 
will be needed to estimate the effect of the lower DM 
density on the formation of low mass DM halos. Finally, 
we illustrate in Fig. [T^] the formation of larger halo NFW 
concentrations in RcBEL, which is in the direction sug- 



gested by recent measurements |24| - |27l \10. 
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Appendix A: Modification of Gadget2 

For the simulation with /3 = we used the original 
Gadget2 public version. For the simulations with /3 = 1, 
we modified the original code as follows: in the PM part 
of the code we added a second loop on top of the original. 
This loop evaluates the density of the DM component 
and transforms it into the scalar potential by convolution 
with the proper green function, given in k space by 



The tree part of the code was modified so that every tree 
node would include a center of charge (scalar charge, 
equals DM mass times /3) in addition to the center of 
mass. A cell-opening criterion was developed, such that 
the relative error in the sum of the ReBEL and gravita- 
tional forces is consistent with the original relative error 
of the gravitational force. The ReBEL force is added in 
a similar way to the gravitational one, that is, by consid- 
ering the first moment. 

To constrain the scalar force relative error we consider 
the worst-case scenario, this is when we calculate the 
force on a DM test particle, located co-linearly with two 
DM particles at the opposite corners of a cubical cell of 
width L. The relative error is of course independent of 
the global interaction strength /? but not on the scale 
length r^. To first non- vanishing order in L/r, where 
r is the distance of the test particle from the center of 
mass/charge of the cubic cell, and up to a factor of 4/3 
the relative error is given by 



(A2) 



which goes back to the usual gravitational relative error 
Cg — L^/r^ for ^ r. Since the particle is given a scalar 
force and a gravitational one, the total relative error is 
constrained by 



et 



CgFg -\- CgFg 



which then gives the new opening criterion as 



(A3) 



(A4) 



Gik) 



fc2 



(Al) 



where e", F^, and F" are evaluated at the nth time-step 
and a ^ 1 is the bounding relative error per time step. 



Appendix B: Simulation tests 

The modified code was checked at the large, cosmo- 
logical scales, by comparing simulations to the ReBEL 
linear theory with a baryonic component included (as de- 
scribed in App. [C]). This is illustrated in Fig. [181 where 
the linear approximation is calculated for the baryonic 
component PS at redshift Z = 100 and Z = Z, and for 
the DM component only at the latter. These calcula- 
tions use no more than the measured PS of the DM at 
Z = 100, a time when the decreasing solutions have al- 
most vanished. The linear approximation seems to fit 
very well the simulations at the linear scales. At smaller 
scales the modified code was checked by two more tests, 
as described in 14^ . The /3 = simulation was checked at 
all scales by comparing the PS shown in Fig. [3]to the PS 
of the corresponding GIF simulation in |l09l |. 
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FIG. 18: Simulation PS compared with linear theory, as de- 
scribed in appendix [Cl for 13 = 1. 



Appendix C: Linear Theory 

First we write the fundamental equations of fluid mo- 
tion in physical coordinates. These consist of the Euler 
equations, 



du. 



VP 



du,. 



dt 



V0,„ + V0,„ = (Cl) 
"^^..r + V0,„ + = 0, (C2) 



where u-^_^^ and u^^^ are the baryonic and DM velocity 
fields, and d/dt is the usual convective derivative, and is 
equal to d/dt + u ■ \7. The symbols (/)^_^_,, 0^^^, and (j)^, 
represent the baryonic, DM, and scalar potential. These 
are given by the poison equations 

- 4nGp,^^ = (C3) 
V2 0^.^, - 47rGp,„ = (C4) 
V'(03-r-')-4^G/3p,„ = (C5) 

where the last equation is a screened poison equation 
describing the scalar potential, and is a function of 
time, and is proportional to the scale factor a{t). Since 
this study deals with scales much larger than the Jeans 
scales, and the simulations are dissipation-less, we neglect 
the baryonic pressure and radiative cooling terms. Thus, 
we skip the baryonic equations of energy and state, which 
play no role in the structure formation dynamics, and 
conclude this presentation with the continuity equations: 



dt 
dt 



(C6) 
(C7) 



We define Su^^^^ as the baryonic peculiar velocity and 
f> as the dark matter peculiar gravitational potential. 



In addition, it is assumed that the baryonic mass fracture 
is much smaller than that of DM (hence the gravitational 
potential is given only by the DM) . Under these assump- 
tions, the linearised Euler and Continuity equations, for 
the baryons, in physical (Eulerian) coordinates may be 
written as 111 101 



(C8) 
(C9) 



where S = Sp^^^^/pbar, and Sp^^^^ is the peculiar baryonic 
density. From now on these equations will be developed 
working in comoving coordinates. To achieve this, we 
transform V — ?► V/o and du, — >■ av, , where a is the 

/ bar bar " 

scale factor, to get 



(CIO) 
(Cll) 



where Vt/)^^^^ is given by the correspondent poison equa- 
tion. In A'-space V ik and — > — fc^ so in trans- 
forming these equations to iiT-space we have 



(C12) 
(C13) 



Eq. IC12l is then multiplied by ik and Eq. IC13l is differ- 
entiated relative to t to get 



S^.. = -ik • 



(C14) 
(C15) 

taking into account Eqs. IC13l and lC15l we have 

-L, =fc'0,„v+2ffC- (C16) 

The poison equation of gravity in comoving coordinates 
is 



2fc2 



(C17) 



with fl^^ = Pa„^Pcrita and pcrit = SH^/8itG. Hence, 
Eq. IC16I can be written as 



2Hi ^ -H^n 6, 

bar 2 



(C18) 



and it is useful to define £7^ in terms of it's current value 
ri^P as VL^ = U^gHQH~^a~^ . in a similar way we can 
derive the linear equation for perturbations in the dark- 
matter density, except that this time the total potential 
that a DM particle feels is the gravitational potential 
plus the scalar potential, giving rise to a modified poison 
equation [l| 



/3 



l + (kr,)- 



(C19) 
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and so all in all, the perturbation equation is 



where 



and 



C= 1 



a 



f3 



1 + (kr. 



(C20) 
(C21) 

(C22) 



where the curvature parameter is fip = S^„q + ^,^„ ■ Tak- 
ing z = 1/a— 1 to high numbers, one can see that fl^ 1 
and the density fluctuations behave as in an Einstein- 
DeSitter universe. 

for a constant in an Einstein De-Sitter universe {fl = 
1 and il\ = 0) a oc t^/"^ and H oc t^^ it can be checked 
by inspection 1] that S^^_^^ oc t'P where 



.^^(25- 



24/3 



1 



(C23) 



in contrast to the p = 2/3 factor when no scalar inter- 
actions are present. As for the baryonic fluctuations, 

In this study we generalize the results of [l| for arbi- 
trary fl^^g and r^Ao' while using regular scaling relations. 
For commodity reasons, from now on wherever we refer to 
(5(57^^, C) if either or C equal 1 they will be omitted. 
So for example, the notation of S{rt^^ = IjC*) becomes 
S{C) and S{V,^^ = 1, C = 1) becomes just S. In addition, 
when no species is defined for (5, it may be either one of 
them (either Sdm or Sbar)- 

We claim by numerical integration, that a good ap- 
proximation for the solution of Eqs. IC18l and lC20l is 

This can be seen in Fig. [121 The relative error of this 
approximation is quite small, as can be seen in the lower 
panel of the same figure. 

for the case of C = 1 we have /3 = and there are no 
scalar interactions - thus we may use standar d app roxi- 
mations for the base of the RHS. for example |111 | 



(C25) 



Since S{C) oc and (5 oc a we can see fror dC24l that 
5{^l^, C) oc 5{Vl^^) ''^^ . And so we get an important result 
for velocity fields, namely 



dlog((5(17.„,c)) 



s^dfogmj) 

2 dlog(a) 



(C26) 



and the derivative term of the RHS can be approximated 
in the usual form, for example [sT 



d log(a) 



0.6 



(C27) 



this is almost independent of A, as shown in [ll2|, which 
also gives a better approximation for a fiat universe: 

M)).o;>^(.-io„,a.o,„,),c.s, 

Next we define the baryonic velocity potential in K- 
space as i'k(f>y = Vj^^^ so Eq. IC13l is 5^^.^^ = fc^i^i,. Noting 
that (5^,,^ = fS^^^H we then have 



(C29) 



where /i is the normal factor /, with C — 1. This result is 
of course good also for the DM fiuctuations. The velocity- 
density relation is thus 



l^PhK.H = -*k . V, 



(C30) 



Linear theory is further extended by giving the baryons 
mass. In this case, Eq. IC20I becomes 

+ 2i/^,„ = \h\^1,^C5,^ + n,.J,J. (C31) 

Lets assume that (5,_,^^ — f with constant /. Then 
the last equation becomes 



-f 2HL ^-H^n^ C*S^ 

dm dm ^ dm d 



(C32) 



where C* = ^7bar/(/^dm) + C. The solution to this equa- 
tion was already given before. As for the baryons, Eq. 
IC18I becomes 

+ 2i?L, = IH\^,J,.^^ + n,J,J (C33) 

Then, using the definition of /, this becomes 

S,^ + '^HS,^ = yH\n,J,^/,f + n,J,J (C34) 



which is true if and only if it becomes Eq. IC321 That 
means. 



(C36) 



(C37) 



which gives 
/ = 

2^^d.„ 

and indeed this factor is constant, since fl^ and fl^ 

' dm bar 

have the same dependence on time. 

This theory is compared to numerical simulations, as 
seen in Fig. [THl and as described in App. |B] This com- 
parison indicates on the correctness of both the simula- 
tions and the linear theory. 
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FIG. 19: Approximate analytic solution to linear theory. Up- 
per panel compares the numerical calculation in dashed red 
lines, to the approximation, which corresponds to the solid 
blue lines. The comparison is made both for a flat universe 
and for an universe without dark energy, for different values 
of C, as indicated in the flgure. The bottom panel shows the 
relative error in the approximation. 



